SOME OBSERVATIONS FOR MEAN-FIELD SPIN GLASS MODELS 



SHANNON STARR AND BRIGITTA VERMESI 



Abstract. We obtain bounds to show that the pressure of a two-body, mean-field spin glass is a 
Lipschitz function of the underlying distribution of the random coupling constants, with respect to 
a particular semi-norm. This allows us to re-derive a result of Carmona and Hu, on the universality 
of the SK model, by a different proof, and to generalize this result to the Viana-Bray model. We 
also prove another bound, suitable when the coupling constants are not independent, which is what 
is necessary if one wants to consider "canonical" instead of "grand canonical" versions of the SK 
and Viana-Bray models. Finally, we review Viana-Bray type models, using the language of Levy 
processes, which is natural in this context. 



1. Continuity of pressure with respect to the coupling distribution 
Let us consider a mean-field spin glass Hamiltonian of the form 

N N 

-Hjv(cr) = ^ ^N{^,j)cri(JJ+h^(Ti, 

where the coefficients JAr(j, j) are i.i.d. random variables, chosen from some distribution represented by 
a cumulative distribution function, Fn, and /i S R is nonrandom. The spins themselves are assumed 
to be Ising type spins, so that a — (cti, . . . ,(Tn) lives in fl^ — {+1,-1}^. This general definition 
encompasses the SK model, where Fn is the c.d.f., for a J\f{0, (3"^ /2N) random variable. That is, all 
the JAr(«, j)'s are normal (Gaussian) random variables with mean and variance 1/2N. This definition 
also includes a nice version of a Viana-Bray type model, wherein Fn is the c.d.f. for a Poissonized 
Gaussian random variable. More specifically, for each i and j, we have 

K(ij) 

Xhj) = ' 

k=l 

where K(z,j) is a Poisson random variable, with mean a/N, and all the random variables gk{i,j), for 
k = 1,2, . . . , are i.i.d. A/'(0, /3^) random variables. 

Note that we write JAr(«, j) and H]s[{(j) using the sans serif font to indicate that they are random 
variables. As usual, we also define the random partition function, 



and the "pressure" 



Pn = ^ln(Zjv). 



This random variable is not actually the pressure, which is only defined in the thermodynamic limit. 
But we will call it the "random pressure" , anyway. Let us write 

which is the definition of the "quenched pressure" . Note that at this point we explicitly denote the 
underlying distribution Fjv for the coupling constants (JAr(z,j) : 1 < i,j < N). This is because 
Pn{Pn) is not random, but a function of Fjy. 

A very basic, but important bound is the following. 
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Lemma 1.1. Let wi,W2, ■ ■ ■ ,wk > be given, and suppose wi + ■ ■ ■ + wk > 0. Also, suppose 
si, . . . ,sk G [—1,1] are given. Define a function / : R — > M by 



fix) = In (^E"'^^""') 



Then f is globally Lipschitz, with Lipschitz constant 1. 
Proof: Note that 



fix) 



OkSk , where dk 



WkC-^ 



k=l 



for k = l,...,K. 



So /'(x) is in the convex hull of {si, . . . , sk}, which is a subset of [—1,1] by our assumption. So 
I /'(a;) I < 1 for all x, which proves the claim. □ 

By the lemma, we see that pjv is a jointly globally Lipschitz function of the random coupling 
constants (Jjv(«, j) : 1 < i,j < N). Let !Fi denote the set of all c.d.f's F which have finite first 
moment. Then, from the above wc sec that, if G J^i, then pjv also has a finite first moment. 
Therefore, y»jv : ^-"1 ^ M is a well-defined function. Part of our goal is to understand the continuity 
properties of pn. Another goal is to derive useful inequalities for the study of "real" mean-field spin 
glass models, such as the SK and VB models. 

1.1. Continuity of pn- Given Fe.Fi, let us define 



ao{F) 



ln(cosh(a;)) dF{x) . 



and for A; = 1, 2, ... , define 



afc(F) = T / tanh'=(a;)(iF(a;). 
Proposition 1.2. Suppose Fj^^F^ S T\. Then 

00 

\pn{Fn)-Pn{F)n\ < W^Fn) - ak{FN)\ ■ 



k=0 



Proof: Enumerate all A''^ pairs {i,j) in any way, as {ii,ji),---,{iN^,jN^)- Then the two Hamiltonians 
associated to Fn and Fn are 



N 



N 



and 



n=l 

For n = 1, . . . define 



i=l 



n=l 



i=l 



N 

'„{a) = exp I ^^N{ik,jk)(Tik(Tj^ + ^N{ik,3k)(Tik(^3k + h'Y^ai 

k=n+l 



, k=l 



i=l 



Then, by a telescoping sum, 

JV2 



9N-PN = 



N 



(1) 



JV" 



SOME OBSERVATIONS FOR MEAN-FIELD SPIN GLASS MODELS 



3 



Let us define w„(cr) = w„((t)/ J^aefi w„(cr). Then we can rewrite 



dn = In J2 w«(^)eJ"(*"'^")"-"- -In J2 w„(a)eJ"(^-^")'^-'^- , 
Vo-enjv / \o-esiiv / 

since the normalizing muhiphers cancel in the difference of the logarithms. Also, let us define 

Then, since e^"^ — cosh(a;) + CTsinh(a;) for ±l-valued variables cr, 

gJ«(»„j„)-.„a,„ ^ cosh(Jw(z„, j„)) [1 + (T,„(T,„ tanh(Jw(i„, j„))] 



In ( w„(CT)eJ' 



ln[cosh(JAr(z„, j„))] +ln[l + tanh(JAr(i„, j„))(cri„(Tj„)f 



Since | tanh(JAr(i„, j„))(o'i„(Tj„)„| < 1, we can use the Taylor expansion of ln(l + x): 



ln(l + x) = 



k=l 



Therefore, defining 



hoix) = ln(cosh(a;)) and = tanh (x) for k — 1,2,. 



we have 



0fe(JAr(»n,in)) " 0fc(JAr(jn, jn)) 



Importantly, by our assumption, JAr(i„,j„) and JAr(i„,j„) are independent of all iN{ik,jk) aand 
^N{ik,ik) for k ^ n. Therefore, SN{in,jn) and SN{in,in) are both independent of {'Ji„<Jj„)n- Also, 
obviously, E[0fe(J)] = ak{F) for an F-distributed random variable J. 

oo 

E[d„] = ao(i^) - ao{F) - ^^(-l)' E [((a,„a,„)„)^l [afc(F) - . 



fe=i 



Since, |(a'i„cr7„)n| < 1, 



|E[d„]| < 51|afc(i^jv)-afc(i^w) 



This is true for all 1 < n < N'^. Plugging into ([T]), we obtain the result. 



□ 



1.2. Re-derivation of a result of Carmona and Hu. In an important paper, Carmona and Hu 
proved universality for the Sherrington-Kirkpatrick model [5] . Using the proposition, we can re-derive 
part of their results, in fact with slightly weaker hypotheses. Suppose is a cumulative distribution 
function such that 

/oo 
xUF{x) = /3V2, 
-OO 

and 

/oo 
xdF{x) = 0. 
-oo 

In other words, letting J be distributed according to F, we have E[J] ~ and Var(J) — /3^/2. The 
c.d.f. of J/a/ZV, then, is Fjm{x) — F{\/Nx). Let us write pAr(i^Af) just to make expHcit the dependence 
of Pat on the distribution of all the random coupling constants (JAr(i, j) : 1 <i,j < N). 
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Let Fp be the c.di. for an A/'(0,/3^/2) random variable and F^^(x) = Fp{VNx). Then, defining 
P^iP) '■= PnIF^ p), this is the quenched pressure for the SK model at inverse temperature /3. By an 
important result of Guerra and Toninelli, the thermodynamic limit, p^^{0) := limjv— >oo exists. 
Among other important results, Carmona and Hu proved the following result, however with slightly 
different hypotheses: 

Corollary 1.3. Suppose F satisfies ^ and and Fn{x) = F{y/Nx) for all N > 0. Then 

Pn{Fn) P^'^iP) , 
as N oo, almost surely, and in probability. 

This is a non-quantitative version of Carmona and Hu's Theorem 1. They obtained explicit bounds 
on \pn{Fn) — p^^iV(/3)|. On the other hand, to do so, they assumed that F also has a finite third 
moment, which is slightly stronger than our hypotheses. (But without such an assumption, it is 
impossible to get quantitative bounds.) We will re-prove CoroUarv 11.31 in order to demonstrate that 
it can be derived from Proposition 11.21 rather than from Carmona and Hu's alternative, and very 
interesting, approach of "approximate Gaussian integration by parts" . The reason we do this is that 
later we will use Proposition 1 1 .21 to prove an analogous result for the Viana-Bray model. (Incidentally, 
Carmona and Hu's result was extended to quantum spin glasses by Crawford in We have not 
considered whether a similar generalization is possible for Proposition 1 1.21 ) 

We can prove convergence in mean, easily, using the proposition. But then convergence in probability 
follows trivially using an important and well-known result of Pastur and Shcherbina. In an important 
paper, they proved that the spin-spin overlap order parameter of the SK model is not self-averaging 
in a certain regime [7]. They also proved the much simpler, but even more widely cited, result in an 
appendix, that the random pressure is self-averaging. It is that second result which we now state, and 
for completeness prove, in the context of general mean-field spin glasses. 

Lemma 1.4. Let Fn satisfy pl2{Fn) < oo. Let p^r denote the random pressure associated to the 
distribution of the coupling constants with c.d.f. Fjy. Then Var(pAr) < a'^{F]\[) = ^2{Fn) — ^i{Fn)'^ ■ 

Proof of Lemma II. 4t The proof follows the standard martingale method. Again, enumerate all N"^ 
pairs {i,j) as (ii, ji), ■ ■ ■ , {iN^ , Jn^)- Let JF„ = a{j{ik,jk) ■ k < n) he the cr-algebra generated by 
J(*fc; jfc) for k < n. Let !Fq be the trivial tr-algebra. For n — 0, 1, ... , N'^, define M„ = E[pAr | JF„]. 
Note that — Mjya while E[p7v] ~ Mq. Therefore, by a telescoping sum, 

(4) Var(pA.) - ^E[M2 -M2„i]. 

n=l 

In addition to the coupling constants JN{i,j) distributed according to F^, let us define additional 
independent copies of these random variables, called ^]^{i,j) and J^(z, j). Define, for a,(3 £ 

n-l N 
k=l k=n+l i=l 

Also define 

P^.n = ^l^^f E exp(-H-'^(a))] . 

\o-enjv / 

Then 

E[Ml - M^J = iE[(p]^^„ - p^^^J(p;^^„ - p^^J] . 
Finally, we want to bound this. So, define 

k=l /c— n+1 *=1 
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Then, for /3 = 1,2, 



I 1,13 2./3 I _ 
|PAf,n Pn,7i\ — ^ 



In ( 



By Lemma Fl-H the right hand side is bounded above by N ^|J]Y(i„, j„) — ■^'Ni'in, jn)\- Therefore, 



Using this bound for aU 1 < n < 7V^, and plugging in to ([4]) gives the resuh. 

Proof of Corollary II. 3t By the triangle inequality and Proposition [T21 we can bound 

\pn{Fn) - Pl^'m < Sn{Fn) + SNiF*^,0) , 

where 



□ 



Sn{Fn) 



Nao{FN) ■ 



Na2iFN) 



f3' 



k=l 
k^2 



We will show that 5n{Fn) — > 0, as ^ oo. Since F'^ also satisfies ^ and ([3]), this also implies 
Sn{F^^p) — > as —> oo. Note that 

^~Nao{FN)= / —dF{x)-j N hi{cos\i{x)) (IFn {x) 

= £(^_„,„(eo.h(-^))).n.). 

Note that, since ln(cosh(a::)) < a;^/2, the integrand is nonnegative. Similarly, 



^-^tanhM^))dF(.), 



and the integrand is nonnegative because tanh^(a;) < x^ . Also, notice that a2k{FN) > for all 
fc = 0, 1, . . . , because tanh^''(a;) > 0, and also 

V ^1^1^ = -1 ln(l - tanh2(x)) - tanh^Ca;) = ln(cosh(.T)) - ^^^^^LM . 

fc=2 



Therefore, 



0"^ 

NaoiFN) - Y 



+ 



Na2{FN) - ^ 



A^^|a2fe(Fjv)| = / ( x^ - N ta.i-ih^ (x/Vn)) dF{x). 
k=2 -^-"^ 



By the dominated convergence theorem, the last quantity approaches as A^ ^ oo. Note that, since 
Mi(F) = 0, 

pOO 

-NaiiFiv) = N{fii{FN) - ai{FN)) = N {x - ta.nh{x)) (IFn (x) 



N 



oo / rx 



oo \J0 



tanh (y) dy ) (iFjv(x) 



So, 



N\aiiFN)\ < 



= N I tanh^(y) 
tanh^(y) 



\x\>Vni 



sgn(a:) dF{x) dy . 



V 



\x\>VNy 



x^ dF{x) dy . 



SOME OBSERVATIONS FOR MEAN-FIELD SPIN GLASS MODELS 



6 



Since H2{F) < oo, this quantity also approaches as iV ^ cx), by the dominated convergence theorem. 
Finally, 



fc=i fc=i -^""^ 

/•°° ^ |tanh2'=+i(:E)| ,^ ^ ^ 



CX3 



^in n+!;""::S";: )-itanh(.)i)df^.(.) 



1 — I tanh(a; 
(|a;| - |tanh(a;)|)dFAr(a;) . 

Therefore, the same argument as the one just above shows that 

oo 

lim Ny^\a2k+i{FN)\ = 0. 

k=l 

This completes the proof that limAr^oo ^n{Fn) = 0. Therefore, it also shows that limAr_,oo |pAf(^Ar) — 
p^{P)\ = 0. Since hmAr_foo = P^^iP), in order to complete the proof, all we need to show is 
that pn{Fn) — pn{Fn) ^ 0, a,s N oo, in probability. But, by Pastur and Shcherbina's bound, 

E[{pN{FN)-pN{FN)f] < pi2{FN) - 

So, we have the even stronger result: pjv(^jv) — Pn{Fn) in L^. □ 

1.3. An application to the Viana-Bray model. Now, let us now reconsider Proposition II. 21 in the 
context of the Viana-Bray model. Define pJP{a, /?) — pn{F^ ^ ^), where 



^N,a,f3 - e 2^ (^1,^) , 

k=0 

in which ★ is the convolution product. Hence, defining K to be a Poisson-(a/Af) random variable, and 
defining gi, g2, • ■ ■ to be i.i.d., M{0, 1) random variables, F^ a p c.d.f. for the random variable 

K 

J -Eg- 

k=l 

This is the "Poissonized Gaussian" coupling used in one version of the Viana-Bray model. 

Corollary 1.5. Suppose that, for each N > 0, the sequence (/at.o, /a^.i, ■ • ■ ) is a probability mass 
function, such that 



lim Nf^ 1 = a and lim N y kf^ fe = . 

Af— >oo A— >oo ^ — ' 

k=2 

For each N > 0, define 

oo 
k=0 

Then 

lim |pjv(i^A)-Pw^(a,/3)| = 0. 

A^oo 

As a particular application, we could take fN.k = (1 ^ "^)'^fc.o + "^'^fe,i- I-^-, instead of taking K to be 
Poisson, with mean a/N, we could take it to be Bernoulli with the same mean. 
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Proof: Note that 

OO 



^ POO 

J2 MFn) - fN,ia2k{Fl0)\ = / 21n(cosh(x)) [dFN{x) - /jv,i dF^pix)] 

/OO 

k=2 

= -y fc fN,k ■ 

k=2 

Defining ^,{a) = e^"/^(a/iV)'''/fc!, we see that exactly the same property is true of it. But also, by 
assumption and calculation, 



hm 7V^fc/i^,,(a) = lim a(l - e""/^) = = lim iV^fc/Ar,^. 

Therefore, 

CSO OO 

limsupTV^ |a2fe(FAr) - a2fc(F^„^)| = limsup N\fN,i - f*N iia)\^a2kiF* 
fc=o ' ' fc=o 



(3^ 

< — limsupiV|/w,i - /^ i(a)| . 

^ N^oo 



But, of course, by assumption and calculation 



lim iV/^i(a) = lim e""/^a = a = lim NfN,i. 

Since ak{FN) — ak{F^ a p) — ^7 because the Gaussian is symmetric, we then see that 

00 

hm NY,\ak{FN)-a,{F;,^^^p)\ = 0, 

k=0 

and we can apply Proposition II. 21 □ 



Also note that, again, Pastur and Shcherbina's self-averaging bounds prove that pn —Pn{Fn) 0, in 
L'^, as N ~* OO. Indeed, 

OO ^2 ^ 

k=0 k=0 

and the right-hand-side goes to 0, as iV oo, by hypothesis. 



1.3.1. Poisson thinning. Let us briefly address one possible point of confusion. The Hamiltonian we 
wrote for the Viana-Bray model was 

ij — l k—1 i—1 

where, for each {i,j) G {1,...,N} x {1,...,A^} the random variable K(j,j) is a Poisson random 
variable, with mean a/N, such that all the variables {K(i,j) : I < i, j < N} are independent, and 
all the random variables gk{i,j), for k = 1,2, ... , are i.i.d. A/'(0,/3^) random variables, all of which 
are independent, and independent of the K(i, j)'s. This is not literally the model that is written down 
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in some references on the Viana-Bray model. Let us call the model we wrote above the Poissonized 
Viana Bray mode^. The original Viana-Bray model considered by many authors is 



N 
i=l 



K 

fc=i 



where K is a Poisson random variable, with mean aiV, and ii, i2, • ■ • and ji, j2, ■ • ■ are i.i.d., random 
variables, uniformly distributed on the N sites {1, . . . , N}. Also, gi, g2, . . . are i.i.d., A/'(0, 0^) random 
variables, independent of everything else. 

These two versions of the model are statistically equivalent. So all expectations of all functions 
of the Hamiltonians are equal, including the quenched pressure. To see this, we use a well-known 
property of Poisson random variables, which is commonly called "Poisson thinning" . We can construct 
a direct correspondence between the random variables of the first Hamiltonian and the second one. 
For instance, given K, ii, \2, ■ ■ ■ , ji,j2, • ■ • and gi,g2, . • . , do the following: First, let K(i, j) — : 
fc < K , ife = j , jfc = i} for each 1 < i, j < N . Note that 



E 



N 



exp ^ AjjK(i,j) 



= E 



exp 



= E 



N- 



N 

E 



exp \ aN \ 



N N 

TE^^- 

i=i ]=i 



N 

nexp(^[e^--l]) 



which is exactly the joint moment generating function for i.i.d., Poisson random variables with mean 
a/N. Therefore, the random variables K{i,j) have identical joint distribution to K(i,j). Similarly, 
there is a way to construct gfc(«,j)'s from the g^'s, by merely letting gkihj) equal gnt(i,j) where 



'^kihj) is the kth smallest integer n such that (i„,j„) — 
see that the collection of {gk{i,j) ■ 1 ^ hj ^ N , k — 1,2, 
{gk{i,j) ■■ l<i,J<N, k = l,2,...}. 



Using independence, it is trivial to 
. } is equivalent to the collection of 



2. Continuity for non-independent couplings 



A key ingredient in the proof of Proposition 11.21 was the assumption that all the coefficients were 
independent. But for some purposes, one wants to drop that assumption, instead assuming that all 
the coefficients and are defined on a common probability space, and are close in some 

sense. Let us state a bound that works in that case. 



Proposition 2.1. Suppose the following random variables are defined on one probability space: a 
random integer N > 0; random spin sites ii, . . . , ji, . . . ,ji\i € {!,..., TV}; and random couplings 
Ji, . . . , Ji\i and Ji, . . . , Jat, which may be dependent and possibly not identically distributed. Define the 
random Hamiltonians 



The 









H 


N 





N 



and 



E 

N 



N 
i=l 



In 



< — E 

- TV 



■ N 

El 

.n=l 



-'^We thank an anonymous referee for suggesting this name, as well as for raising the issue of demonstrating the fact 
that the two versions of the model are statistically equivalent. 
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Proof: Consider the linear interpolation H^ t = tHj^ + i 
and integrating over (0,1) it readily follows that 



dt 



1 - t)HN. Then 

^ — Hjv,t(o' ) 



(Hjv(a) - Hw(a)) 



1 



N' 



E 









)1 




H 




N 





In J2 ^'""^"^ 
Now conditioned on N, for all a ^ il^ , 

Hw(a) - Hjv(a) 
Taking expectations gives the desired result. 



< 











)1 




< — E 


max 






- N 





Hw(a)-HAr(a) 



n -'n 



□ 



2.1. A canonical, versus grand canonical, version of the Viana-Bray model. An equivalent 
definition of the Hamiltonian for the Viana-Bray model [9] is 



N 



-HAr(tT) = y Jn fJi„ q'j„ + y (Ji , 
n— 1 z— 1 

where: K is a Poisson random variable, with mean aN; the J„'s are i.i.d., A/'(0,/3^) Gaussian random 
variables; and ii, 12, ■ ■ ■ Ji,j2, • ■ ■ are i.i.d., integer-valued random variables, uniformly distributed on 
the set {l,2,...,iV}. Then, 



In J2 '''""^'^ 



Note that the number of edges present is a random variable K. 

In statistical mechanics, we often consider the canonical ensemble to be a model of a gas where the 
number of particles (but not the energy) is held fixed. In the grand canonical ensemble, the number 
of particles is random, although usually highly concentrated around its mean value. Since the number 
of edges in the Viana-Bray model is random (but highly concentrated around its average), this is like 
a grand canonical ensemble for the number of edges. Consider an alternative "canonical ensemble" 
Hamiltonian 

[aN] N 
n— 1 i—1 

where all the J„'s are i.i.d., Af{0, 0^) Gaussian random variables, and ii, i2, . . . , ji, j2, ■ • ■ are as before. 
Now the number of edges is nonrandom: it is \_aN\ is nonrandom, the greatest integer < aN . Define 



In ^"^"^'^ 



Corollary 2.2. With the definitions above, 



Remark 2.1. A qualitative version of this result was stated without proof in the papers of Franz and 
Leone [5] and Guerra and Toninelli [5] . 
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Proof: Define gi,g2, ■ ■ • to be i.i.d. Af{0, 1) random variables. Let K be a V{aN) random variable, 
and let ii, ■ ■ ■ , ji, j2, ... be i.i.d., and uniform on {1, . . . , N}, as before. For n < K, let J„ — /3g„. 
For n > K, define J„ = 0. For n < laN\ define J„ — (3gn- For n > [aiVj define J„ = 0. Let 
N = max(K, [aTVj). It is easy to see that Hjv((j) and Hjvlc), defined as in Prop 12. l1 have the correct 
distributions for the Viana-Bray and the "canonical ensemble" models, respectively. Therefore, 

N 



\pTia,P)-pJP{a,P)\ < 1 



< A(i + Var(K)i/2) Var(g^)i/2 



^(1 



/oN). 



□ 



2.2. The canonical, versus grand canonical, version of the SK model. An analogous result 
holds for the SK model. A definition of the SK model is as follows: Let Xjv be a xn'^ random variable. 
Let V^r = : 1 < i,j < N) be a uniform random point on the unit sphere 



jW-i 



Then, defining Jjv(i, j) = /3Vw(i, j)Xw/V2iV, 



N N 
a^UN — 1 ^—1 



Again, note that the norm of the coupling constant vector, Jat = (Jjv(*,i) : 1 < < N), equals 
Xn/\/2N, which is random itself. One could consider this to be a grand canonical ensemble. In the 
"canonical ensemble" the only thing random about the coupling constant vector would be the direction. 
Therefore, define — /3V(i, j)-\/iV/2, and define 

N N 



This is the analogue of the SK model, but where the couplings constant vector, (Jjv(i, j) : 1 < «, J < 
A''), is constrained to lie on a sphere with radius Rn and satisfying i?^ = (PN/2 = E[|| Jjv||^]. 

Corollary 2.3. With the definitions above, 



1 



Proof: By direct application of Prop 12.1 



1 



\pT{i3)-PTm < 



N 



2J = 1 

AT 



< 



V2iV3 



J2 E[|V^(z,j)|.|XAr-A[ 



N 



V2iV3 
V2N 



2^1/2 



i,i=l 



l[X%] +N^- 2AE[X7v]) 



1/2 
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But, as Xat is XAf2, that means E[X^] = N"^ , while 



Therefore, 

E[X^] + TV^ - 2ArE[Xjv] = 0(1) . 

□ 



3. Mean-field spin glass models with infinitely divisible couplings 

All of the results so far were motivated by the beautiful results for the Viana-Bray models obtained 
by Franz and Leone |5] , Guerra and Toninelli [B] , and De Sanctis [3] . Additionally, we were motivated 
by Carmona and Hu's universality result for the SK model [2]- The Viana-Bray model, as studied in 
the papers listed above, basically relies upon one property of the coupling distribution. That is that 
Fm is infinitely divisible. Let us digress briefly, to discuss infinitely divisible distributions. 

We will specialize our attention to symmetric distributions. Suppose that A is a nonnegative measure 
on (0,oo) (not including 0) satisfying 



oo 



mm 



{y\l)K{dy) < oo . 



0^ 



(We write in place of J(^q ^y) Also suppose w is a nonnegative number. Then one can define a 
function 

(K.v){k) = — + {1 - cos{ky)) A{dy) . 

This function is conditionally negative semidefinite, and at 0. In other words (c.f. Schoenberg's 
theorem), exp [— ^(^^ ,;)(fc)] is a positive semidefinite function of fc G R, and equals 1 at fc = 0. 
Therefore, by Bochner's theorem, it is the characteristic function of a unique c.d.f. We define i^(A,t)) 
to be this c.d.f. Thus, 

oo 



(A,«)l 

-oo 

It is a basic fact that F(Ai,i>i) *-P"(A2,ti9) = -P'(Ai-i-A2,tii+i'2)- Therefore, Fi^i^^^-^ is infinitely divisible: in 
fact i^(A,i,) = P[\/n v/n)- Levy-Khiuchinc formula, specialized to symmetric distributions, every 

symmetric, infinitely divisible c.d.f. is of this form for a unique pair (A, v). 
Let us suppose that we have a spin glass Hamiltonian, defined as previously 

N N 

-Hjv(ct) = ^ Jjv(», j)q-iO-3 + fey^gj ■ 

— l 2—1 

But, now we suppose that the random couphngs, (JAr(i,j) : 1 < i,j < N) are i.i.d., and distributed 
according to F^ = F{A/N,v/N). We may denote F := Fi = F{A,v). For a pure Gaussian, with 
V = /3^/2, we have Fjy = F{0, /3'^ /2N) = -F^^, as before. For the Poissonized Gaussian, we have 

Pn,c.,0 = Fn= F{aA*plN, 0), where 

Ki{dy) = — ^=^i(o,oo)(y)d?/. 



Therefore, this does, indeed, generalize the two cases we considered before, of the SK model and one 
version of the Viana-Bray model. Then we write 



In ^"^"^"^ 
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We will not necessarily introduce a new symbol for the Hamiltonian, when the underlying distribution 
for the couplings, : 1 < i,j < N), changes. Rather, we will endeavor to write the distribution 

explicitly, when we take expectations, as in Let us write E(A_t,)[-] instead of E-'^'A'") [•]. 

Let us introduce the fundamental definitions of Franz and Leone, Guerra and Toninelli, and De 
Sanctis. First of all, given a function of n spin configurations, u : (Qn)"" — > R, let us write 

(u) = . . . , (Z„)-i J2 e-[H"('^''')+-+^«('^'"')lw(a(i), . . . , a^")) . 

o-(i),...,(T(")enjv 

Note that, in our choice of convention, this is still a random variable depending on the underlying 
coupling constants. But E(A,t,) . ■ . jCr^"^))] has had the expectation taken (with respect to the 

i.i.d. product of F(^a,v) distributions). Let us also define the degree-n multi-overlap function RN,n '■ 
{QnT ^ M, as ' 



1 " 



N 

(1) (") 

N 

1=1 



Then a fundamental result of Franz and Leone and Guerra and Toninelli, suitably generalized to the 
present context, gives an integral formula for pJ^(A2, V2) — p^(Ai, in terms of the expectations of 
these multi-overlaps. Let us generalize the definition of a2k{F) as follows: 

V f°° 

a*o{A,v) -.^ -+ ln(cosh(y))A(d?/); 



2 J0+ 

V 1 



2 2 



al{A,v) - + - / tanh (y)A(dy); and 



0+ 



a*2k{A,v) = a;^{A) -.^ ^ I ts^nh'^y) A{dy) , for fc = 2,3,. 



2k 



0+ 

1 A" 



(In some sense, when we write (A, v) this really means the distribution A + ^^q, acting on smooth test 
functions cj) satisfying (j){Q) = 0'(O) = 0.) Then the result is as follows. 

Proposition 3.1. Suppose that (Ai, vi) and (A2, V2) are parameters from the Levy-Khinchine formula, 
satisfying the further requirement that ao(Ai, vi) and ao(A2, V2) are finite. Then 

PAr(A2,W2) -Pw(Ai,wi) = aQ(A2,W2) - a5(Ai,i;i) 

°° /"i r / 2\ " 

- X]["2fc(A2,W2) - a2fc(Ai,'yi)] J E(4A2-(-(l-t)Ai,to2 + (l-t)«i) (^(i?Ar,n(CT*^\.--,Cr("))) ^ dt . 

Remark 3.2. One clearly sees the motivation for Provosition \1.2\ in this formula. 

An immediate corollary, along the lines of Proposition 11.21 can be deduced from this. Namely, 

00 

\p*NiA2,V2) ~ PNiAi,vi)\ < ^ |a2fc(A2,'y2) - a2fe(Ai,fi)| . 

This gives a bound which is uniform in N , and is often easier to calculate. Guerra and Toninelli were 
the first to write this type of bound, when they used it to give a very simple proof that 



lim sup 



pl^(a,-^]-p%HP) 



0, 

for all /? > 0. This is known as the "infinite connectivity limit". 

The main advantage of Prop [XT] is that it is an exact formula. For example, Franz and Leone used 
a close analogue of this formula to prove that, for the Viana-Bray model, the thermodynamic limit 
of the pressure exists. Guerra and Toninelli used similar methods to control the high-temperature 
and low-connectivity regions of phase space, demonstrating replica symmetry in that domain. And 
De Sanctis used that method and other arguments to prove an extended variational principle, thereby 
generalizing the results of Aizenman, Sims, and an author [T], from the SK to the Viana-Bray model. 
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A specialized version of Prop 13.11 applicable to the standard Viana-Bray model, is contained im- 
plicitly or explicitly in each of the papers J5] , [S] and [?] . Accordingly, the reader may find the relevant 
proofs, there. However, a new issue arises in the generalized context. Namely, we should prove that 
the pressure function, Pat, is still in the domain of the "generator", despite the fact that it is not a 
typical test-function (because it does not vanish at infinity). Next, we will present the definition of 
the generator, as well as this technical result. 

3.1. The generator. An important fact is that one can define a Levy process associated to the 
infinitely divisible distribution -F'(A,t))- This is a stochastic process (Xj : t > 0). Among other properties 
are these two: for each s,i > 0, the increment (}^s+t — Xg) is independent of JF,, where Jv, is the a- 
algebra generated by (X^ : < r < s); and {Xg+t — Xg) has the c.d.f. i^(tA,tu)- In particular, the 
increments are independent and stationary. (There are also continuity properties of the Levy process 
which we will not need.) See [8], for example, for a reference. Let Cq(U.) denote the set of function 
/ S M which are twice continuously differentiable, and such that 

lim f{x) = lim f{x) = lim /"(x) 0. 

|a:| — *oo \x\ — ^oo |a;| — >oo 

Then the following is a specialization of Theorem 31.5 in 'W/. If / G Cq(R), then E[/(Xf)] is differen- 
tiable, and 

(5) ^'^[•^(^*)] = IE[©(A..)/(X*)] , 
where ©(a.d) is the generator 

^iA..v)f{^) = If'ix) + 1^ + y)- fix) + ^fix - y)j Aidy) . 

In particular, if we instead consider X to be distributed by -F'(A,ti) j ^nd denote the expectation with 
respect to F(^a,v) as E(a,«) , then we have 

^E(tA,t..)[/(X)] = E(tA,t.)[6(A,.)/(X)]. 

The technical fact we want to prove now is the following: 

Lemma 3.3. Suppose (Ai,i;i) and (A2,ut) satisfy ag(Ai,Di) < oo and a2(A2,U2) < oo. Also suppose 
that f : M. ^ R is a function in C^(]R) such that f and f" are in L°°(M). (In particular, then, f is 
globally Lipschitz.) Then ^(Ai,vi}f '^i^d ®(A2,t>2)/ '^'^^ both well-defined and 

(6) E(A2,t,2)[/(X)] - E(Ai.i,i)[/(X)] = / E(4A2+(i-t)Ai,ti,2+(i-t)^'i) [(®(A2,«2) - ®(Ai,t,i)) /(X)] . 

Jo 

This lemma proves that the generator is applicable even to the function pat, which is a Lipschitz 
function of the random coupling constants J(i, j). The beautiful calculation of the generator for this 
function can be read off of any of the references mentioned before, [51 [SI ?] . We will spend the rest of 
this section proving Lemma 13.31 



3.2. Proof. First, supposing aQ{A,v) < oo, let us show that F^^v) finite first moment. This 
is equivalent to checking that ln(cosh(a;)) is integrable, because ln(cosh(a:)) ^ |a:| as |x| oo. But 
ln(cosh(a;)) — e"-'^ ln(cosh(ea:)) is a nonnegative function, for every < e < 1, and it is in Cq(M). So 
the generator apphes to it. But it is easy to see that, defining Ue{x) = e^^ ln(cosh(ea;)), 

/ N 1 , ( sinh^(e-!/)\ . , 

There are several things to note. First, the formula is well-defined. Second, an upper bound is obtained 
by bounding cosh^(ea;) > 1 and (2e)^^ln(l + sinh^(e?/)) = e^^ ln(cosh(e?/)) < ln(cosh(j/)). Therefore, 
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^{A,v)Ue{x) < ao(A, ii), for all x. But finally, by the DCT, linie^o ®(A,i;)We(a;) — > 0, for all x. So, by 
various trivial applications of the DCT, 

IE(A,t;) [ln(cosh(X))] = limE(A_„) - Ue)(X)] = Ihn / ^(tA.tv) [®(A.i,) ("i - '"£)(Xt)] dt 



E(tA,to) [©(A,i;)Ul(Xt)] dt 

A quantitative version shows that E(a,i>) < (1 + aQ(A, v))e. In particular, by assumption ao(Ai, vi) 
and a5(A2,i'2) are finite, therefore, Fi^tAi+{i-t)A2.tvi+(i-t)v-,) has a finite first moment for all t G [0, 1], 
and moreover it is bounded by (1 + max{ao(Ai, wi), ag(A2, W2)})e- Since / is globally Lipschitz, this 
means it is integrable against F(^tAi+{i-t)A2,tvi+{i-t)v2)- 

Now we simply introduce a smooth cut-off, and perform basic estimates. Suppose that ip{x) is 
any function which is twice continuously differentiable, compactly supported, and such that ip{0) = 1. 
Then, defining, fe{x) ~ f{x)ip{ex) this is in with compact support. So equation ([6]) holds with / 
replaced by f^- Since /e — > /, pointwise, and since |/e(a;)| < HV^Hool/Ca^)!, we have 

E(A„..)[/(X)]-E(A„„,)[/(X)] = lim(E(A„.,)[/.(X)]-E(A„„,)[/,(X)]) . 

To prove ^ for / (instead of /e), we just have to show that the right hand side of this equation 
converges to the right hand side of ([6]) as e — > 0. Note that 

\Mx)-My)\ < (||.nioo + e||^'||oo)|x-y|. 

So, for some A < oo, 

lMx + y) + ^Mx~y)- Mx) < A\y\, 
which we use for \y\ > 1. Also, 

< liriloollV-lloo + 2e||/'||^,||^'||oo + e'||^"||oo|/(x)| . 
So, by Taylor's theorem, for some B,C < oo, 



\h{x + y) + \Ux-y)-h{x) 



< {B + eC\f{x)\)y\ 



which we use for \y\ < 1. So, since k\n{cosh{y)) < mm{y^, \y\) < iirin(cosh(2;)), for some < k < K < 
oo, we have 

lfeix + y) + lUx-y)-Ux) 



< (B-FeC|/(x)|)ln(cosh(y)), 



for some different constants B,C < oo, and all j/ e K. Since aQ{Ai,vi) and aQ(A2,U2) are both finite, 
this upper bound is integrable against Ai and A2. Therefore, for some other constants A,B< 00, 

(7) \e^A^^,,)Ux)\ < A + B\f{x)\ for 1 = 1,2. 

But also note that. 



lim Ifeix + y) + l:Mx - y) - fe{x) 

e^O 2 2 



lf{^ + y) + lf{^-y)-fi^): 



for each a; g R and y > 0. Similarly, limj^o f'Jix) — f"{x), for each x e M. So, by the DCT, we know 

l"'^®(Ai-''.)/<^(^) = ^{A,,Vi)f{x) , 

for 1 = 1,2 and all x e M. But also, the upper bound of ([7]) is integrable against i^(fAi-i-(i-t)A2.tfi+(i-t)t)9)j 
for all t, and is independent of e. So, by the DCT again. 



lim ] 



-(tA2-(-(l-i)Ai,it)2 + (l-t) 



vi) [(®(A2,i>2) - ®(Ai,t,i)) /e(X)] 

= IE(fA2 + (l-i)Ai,ti'2 + (l-t)«i) [(®(A2,D2) ~ ^{Auvi)) /(X)] 
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Of course, the integrated version of this is also true, once again by DCT. 

4. Conclusion 

We considered two different types of bounds for the difference of two pressures of two spin glasses. 
In Section 1, we considered a bound which demonstrates that the pressure is Lipschitz with respect 
to a seminorm, ||F|| = J2^=o W2k{F)\. This was strongly motivated by Carmona and Hu's proof of 
universality for the SK model, but generalized to also apply to the Viana-Bray model. In Section 2, we 
considered a bound which is useful if one does not assume that the coupling constants are independent. 
This proved a different type of universality, which was noted, but not proved, in papers of Franz and 
Leone, and Guerra and Toninelli. In Section 3, we briefly reviewed the theory of infinitely divisible 
distributions, and applied it to the Viana-Bray model. All the results in this letter are simple. But we 
hope they add something to the growing wealth of knowledge for mean-field spin glass models. 
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